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We give an exact solution of the bD Einstein equations which in 4L' can be interpreted 
as a sphericahy symmetric dissipative distribution of matter, with heat flux, whose effective 
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' and temperature are related by the Stefan-Boltzmann law. The solution admits a homothetic 
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1 Introduction 



In recent years, there has been an explosion of interest in solutions of the five- dimensional field 
equations of general relativity. This is due mainly to two things. Physically, it is now gener- 
ally acknowledged that the unification of gravity with the interactions of particle physics is best 
approached through the extension of spacetime to higher dimensions, of which (noncompactified) 
Kaluza-Klein space with five dimensions is the basic example. Mathematically, it is now accepted 
that Campbell's theorem guarantees that any solution of the four-dimensional Einstein field equa- 
tions with matter can be derived from a solution of the five-dimensional Kaluza-Klein equations in 
apparent vacuum, provided an appropriate identification is made for the energy-momentum tensor 

Hl-ISl- 

In the present paper, we wish to focus on one situation which is fundamental to astrophysics 
but has not received much attention. Namely, that of a non-static ball of radiation-like matter 
which is spherically symmetric in ordinary 3D space, and whose properties are described by an 
exact solution which has self-similar symmetry. 

The latter symmetry has been extensively applied in studies of fluid flows in accordance with 
the equations of Newton and Einstein |E]-|HI- It's usefulness lies in the fact that when the equations 
are algebraically complicated, it is sometimes possible to find solutions which depend on a ratio of 
the independent variables (rather than on these separately), and that such solutions are frequently 
realized to a good approximation in the real world. This correspondence arises, for example, when 
a solution does not contain any arbitrary length scales, because the physical situation it describes 
does not contain a sharp boundary P]-[T^. 

While the self-similar symmetry has been extensively applied to obtain exact solutions of the 
4D Einstein equations of general relativity, it has not been much studied in the context of the 
5D Kaluza-Klein equations. These are currently the topic of much research, under the names of 
induced-matter and membrane theory. In the former approach, the fifth dimension is responsible for 
the matter/energy content of the universe whereas in the latter approach, the fifth dimension 
provides a singular surface around which are concentrated the particle interactions of the universe 
|,15j . It has recently become clear that the field equations of these two approaches are basically the 
same [12], so in the present paper we will adopt a generic approach, and elucidate the physics of a 
self-similar exact solution of the 5D field equations. 

These field equations are normally taken to be those obtained by setting the higher-dimensional 
Ricci tensor to zero. In this regard, they are the extensions of the Einstein equations of general 
relativity, as verified by the classic solar-system tests of that theory. However, the 5D equations 
are richer than the 4D ones, which means that Birkhoff's theorem does not carry over in simple 
form to 5D from 4D ^7]. Nevertheless, the 5D case is known to agree with the classic and other 
tests of astrophysics, and also gives models which are in agreement with data from cosmology jl8]. 
This concordance, along with the promise of new physics, underlies the current upsurge in interest 
in exact solutions of the 5D field equations. 

However, a gap in our knowledge exists as regards the subject of this paper, namely the case of 
a non-static distribution of matter which is spherically symmetric in 3D space and has an equation 
of state of the hot or radiation type. Such a solution can be expected to be relevant to the early 
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universe, especially in regard to the growth of perturbations in the early stages of inflation and to 
the related topic of galaxy formation. 

Other solutions of this type exist in the literature They are usually derived by solving 

the 5D field equation Rab = 0, and then using the standard technique in which the solution is re- 
expressed as one with Gais = T^p, where the last object describes the effective properties of matter 
[in] in accordance with Campbell's theorem (here Rab is the 5D Ricci tensor, and Gq,/?, T^p are the 
4D Einstein and energy- momentum tensors). The details of this procedure are by now well known 
to researchers in the field, and algebraically tedious. 

In the present paper we give a solution of the 5Z) equations which in AD represents an expanding 
spherical distribution of matter whose effective equation of state is that of radiation. The solution 
may be verified by a fast software package such as GRTensor [20], and the long form of the field 
equations may be found in the literature jEj, [21]. 

Our solution has some unique properties: (i) it is nonstatic and spatially nonuniform; (ii) it is 
curved in 5D, but (iii) admits a homothetic killing vector, and (iv) exists for both signs of the extra 
part of the metric. In these regards, it is different and more realistic than other solutions which 
have appeared, and may in the future deserve more intensive study. 

2 A solution for an expanding ball of radiation 

For an early universe, we are interested in solutions that describe expanding fluids with an equation 
of state which is radiation-like. These should be spherically symmetric in ordinary 3D space. 
Therefore, we choose the line element as 

dS'^ = e'^df - e^dr^ - R\de^ + sin^ed<f) + e^^dy^ (1) 

where y is the coordinate along the extra dimension and t, r, 9 and (j) are the usual coordinates for 
a spacetime with spherically symmetric spatial sections. The factor e in front of $^ can be either 
— 1 or +1 depending on whether the extra dimension is spacelike or timelike, respectively. In what 
follows we choose units such that the speed of light and the gravitational constant are c = 1 and 
SttG = 1. 

In (0) the metric coefficients A, R and $ are, in principle, functions of the extra coordinate. 
However, we assume from the beginning that i/(r, t), A(r, t), R{r,t) and $(r, t) are not functions of 
y. In this way we assure the radiation-like nature of the effective matter in 4D 

The metric in 5D is a solution of the field equations |T3] , which in terms of the Ricci tensor read 
Rab = {A, B = 0,1, 2, 3, 4). Here, for brevity we omit the details of the derivation and state the 
solution 

± C2r2("+i)t2("+3)/"rf^2^ ^2) 

which may be verified by computer. Here a is a dimensionless parameter in the range < < 3, 
whereas B and C are arbitrary constants with the dimensions and L~^" +2a+3)/a^ respectively. 



dS^ = B' 



3r^ 



2^2 



dr - rdr^ 



(3 



■{dO'^ + sin' 
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The parameter a is related to the time evolution of the extra dimension. Indeed, from Q it follows 
that 

^ dt at ' ^ ' 

Since < < 3, it follows that either — <a;<OorO<Q;< \/3. In both cases (a + 3) is 
positive. Therefore, for a negative (positive) a the size of the extra dimension decreases (increases) 
monotonically with time. 

For the four-dimensional interpretation of the solution we identify our spacetime with a hy- 
persurface, located at some value of orthogonal to the extra dimension. This is the common 
assumption in STM and brane theory. Despite the differences between them, the effective matter 
content is the same in both approaches PH]- The appropriate formulae for the calculation of the 
energy-momentum-tensor T^^, for the metric (0) can be found in Ref. [13]. Using ((21) we get 



Tl 



2a2 



2 ^ (3 - a') 3 
2 35 W 3, 

2 

We emphasize that Campbell's theorem guarantees the validity of Einstein's equations G^^, = T^^. 
For the physical interpretation of the spherical distribution of matter Q we notice that: (i) for 
7^ 1 the principal stresses are unequal, (ii) there is a radial flux of energy Tqi ^ 0, (iii) it satisfies 
the energy conditions 

T°>0, (a2<3), T^>\Tll (^ = 1,2,3), (5) 
and (iv) the trace of the energy-momentum-tensor vanishes, viz., 

r° + Tl + r| + T| = 0. (6) 

This is the usual relation for ultra-relativistic particles or photons, and in terms of the average 
pressure p = —(Tl + T^ +Tl)l'i and density p = Tq is commonly written as p = 3p. 



3 Model for the matter 

In order to ascribe a precise meaning to the components ^ it is necessary to adopt a model for 
the matter. The general form of the energy-momentum tensor for anisotropic matter with energy 
flux can be covariantly given as 

T^,u = (p + P±)Uf,u„ - p^g^^ + (p - PA.)XtiXv + + g^M^, (7) 
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where is the four- velocity, is a unit spacelike vector orthogonal to m'^, p is the energy density, 
p is the pressure in the direction of x^t, and pj_ is the pressure in the two-space orthogonal to X/j,- 
The energy flux is described by the 4-vector and obeys g^u^ = 0. 

We now proceed to discuss the matter quantities entering (0) in the comoving and non-comoving 
frames. 



3.1 Comoving frame 

For a spherical distribution of matter, in a comoving frame we have 

u, = 6y/^ x, = 5y'\ q, = 5y'\ (s) 

where q = {—q^iq^Y^'^ ■ Consequently, for the metric (j2I) the matter quantities are as follows 

_ 2 _ 2a^ _ (3 - a^) 

and 

2a 

We see that for a"^ = 1, the radial and tangential stresses are identical to each other p = i.e., 
it is a perfect fluid with an equation of state p = 3p. For any other 0^2 7^ 1 in the range ^2 < 3 
the density and principal pressures are positive, p = p + 2p±. Besides, p > \q\, which is physically 
reasonable. In the present case, the expansion B = is given by 

We note that (a/B) is a positive quantity because the four velocity is future-oriented, i.e., Uq and u° 
are chosen to be positive. Consequently, the fluid under consideration is in continuous expansion. 



3.2 Non-comoving frame 

Thus, in the comoving frame the difference — 1| can be used to measure the "degree" of anisotropy 
of the fluid. We now proceed to show that for any value of a, in the allowed range, one can always 
"jump" to a moving frame where the fluid shows isotropic pressure. 

For a spherical distribution of matter, in a radially moving frame the four-velocity and energy 
flux can be written as 



= (e-'^/^coshtu, e^^/^sinhtu, 0, 0), 



= q {-e^/^sinhuj, e^/^coshoj, 0, 0), (12) 
where is a parameter that measures the radial velocity v, viz.. 



1 V 
coshiu; = —^===, sinhu; = ^===. (13) 



V 
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The introduction of the hyperbohc functions assures the fulfillment of the conditions m^m'^ = 1 and 
M^g^ = 0. 

For perfect fluid {p = p±), the energy-momentum tensor ((Zj) with components in the frame 
yields 



P 
P 



{T^ + Tl-Ti 



(3 - g^) 



P 
3 



(14) 



and 



q = Cp, (15) 

where for simplicity we have introduced the dimensionless quantity (. The explicit expressions for 
Tq , T/ and Tq from ((7j) are 

(9 - a^) 



4 cosh^ — 6C sinh u cosh ci; 
4 sinh^ CO — 6C sinh cosh u 
3( (cosh^ u! + sinh^ a;) — 4 sinh a; cosh u) 



(3-a2)' 

3(a^ - 1) 
(3 - a2) ' 

(3-a2)' 



(16) 



respectively. These equations can be regarded as a system of equations in the regions < a < 1 
and — -\/3 < a < — 1. The solution is 



\a\ 



a 



73' 



coshtu 



9 + — 6|a| 



Consequently, we find 



\ 2(3 

|a|(3-a2) 



sinhtu 



3 + 3a^ — 6|a| 
\ 2(3 -a2) • 



(17) 



3 + 3a^ — 6|a| 



\ 9 + a2-6| 



a\ 



Thus, for a = ±1 the radial velocity vanishes and we recover the expressions (jH)) and ()10|1 for perfect 
fluid in the comoving frame, as expected. It should be noted that the matter quantities are not 
affected by the signature of the extra dimension or the sign of a. The latter affects the evolution 
of the extra dimension Q. 

3.3 Thermodynamic relations 

In order to ensure that a model with heat conduction is physically acceptable, besides the usual 
energy conditions, an appropriate set of thermodynamic relations must be satisfied. These are: (i) 
baryon conservation 



I -,11 



0, 



(19) 
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where n is the particle density; (ii) the first law of thermodynamics 



Td{S/n) = d{p/n) +pd{l/n), (20) 
where T is the temperature and S is the entropy density; (iii) the temperature gradient law 

9m = - uX)Tu + Tu,,uu% (21) 

where k is the thermal conductivity which must satisfy 

> 0, (22) 

and (iv) positive entropy production 

S^, > 0, (23) 

where = Su'' + qf'/T. 

The model under consideration here satisfies all the above conditions. Indeed, integrating ()19|). 
with from (fT^ . we obtain 

n = -4^, (24) 

where no is a constant. We see that the number of particles decreases with time, at any given 
location, which is consistent with the expansion of the fluid. 

The temperature T occurs in (^01) as an integrating factor. In our case, for ()14|) and ()24|) we 
find this factor to be T ~ {rtY ■, where j3 is some constant. If we take (3 = —1/2, then the density 
and temperature are related by the Stefan-Boltzmann law. Namely, 

P = B\3 - a') (^^)\ (25) 

with 

T = (26) 



where To is a constant with the appropriate units. We note that [Tq/B) must be positive in order 
to ensure T > 0. In addition, in view of the fact that the temperature is a decreasing function of 
the radial coordinate, the heat flux is directed outward. 
For the thermal conductivity k we find 

^ " 3\a\V3H{l + av/V3) is) {j^l ' ^'^^^ 

Since (B/Tq) and {a/B) are positive quantities, it follows that k is positive in the allowed range of 
a, which is consistent with ()22|) . 

At this point we should mention that condition ()23j) is ensured by conditions ^2\\ and ()22|1 so it 
is not a separate condition. Therefore, the conclusion from the above discussion is that our model 
satisfies appropriate physical requirements in the regions < a < 1 and —a/3 < a < — 1. 
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4 Properties of the solution 



The solution under consideration has a number of interesting mathematical properties. 
Firstly, the line element (jSJ admits a homothetic Killing vector r]^ in 5D. Namely^, 

^rjQAB = '^QAB-, (28) 

where 

r/^ = -(t, r, 0, 0, y). (29) 

Z a 

Consequently, the spacetime part of inherits the so-called self-similar symmetry, where all di- 
mensionless quantities in the solution can be put as functions of a single independent variable 
jHl-[I3I- Specifically, 

A^Ap = 2^Ap, (30) 

where 

e = lit, r, 0, 0). (31) 

The self-similar nature of the spacetime metric becomes evident in the canonical coordinates t, f 
given by 

t = VAB-H^^-^\ r = VAB^''-^^f\ b = ^=L=, (32) 



a 



where A is some dimensionless constant. Indeed, in these coordinates the 4D metric becomes 

a'^ (3 — a^) 

where the similarity variable ^ is 

? - (3 . (34) 

It should be noted that the five- dimensional solution is not self-similar (in the sense mentioned 
above), only the AD part is. Self-similar symmetry in general relativity generalizes the classical 
notion of similarity, which has successfully been applied in classical hydrodynamics to the description 
of shock waves created in explosions [HI, |Z1- The physical relevance of self-similar solutions is that 
they are singled out from a complicated set of initial conditions [8J. 

Secondly, the present solution differs from others in the literature in several ways: 

1. Solution Q exists for both signs of the extra part of the metric, whereas most solutions that 
have drawn attention only satisfy Rab = for one sign or the other. For example, the standard 
5D cosmology due to Ponce de Leon ^2] only exists for signature (+, — , — , — , — ), while the 
wave solution due to Billyard and Wesson j2SI only exists for signature'^ (+, — , — , — , -|-). 



^Here we use C^gAB = gAB,cV^ + 9cbv'^a + 9acV% 

■^For a recent discussion of the signature of the extra dimension see 
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2. Many of the solutions in the hterature are curved in AD but Riemann-flat in 5D (see the 
catalog of Ref. [6], pp. 95 and 109). However, metric (0) is not flat in 5D. Indeed, Rabcd 
has 12 nonzero components. 

3. The signs of the spacetime components of the Riemann tensor are the same for both signatures 
in (0), whereas the signs of the components associated with the extra dimension reverse 
depending on the sign choice in (|^. 

4. The 5D-Kretschmann scalar K = RabcdR'^^'"^ is unaffected by the signature. For Q we 
find K ~ so it is singular on the hypersurface rt = 0. This kind of singularity is 
common is self-similar solutions. Physically this is because such solutions are applicable far 
from the origin and boundaries, where the initial conditions are unimportant |Hj, [HI, [Zj- 

5 Summary and conclusions 

The line element ((21) is an exact solution of the five- dimensional Einstein field equations in an empty 
5D space. By Campbell's theorem it can be interpreted as a solution of the AD field equations with 
energy-momentum Q. Since the metric (j21) is independent of the extra coordinate, the trace of the 
energy-momentum tensor is zero and the matter distribution is radiation-like. 

In the comoving frame the source can be modeled as an anisotropic fluid with heat flow in the 
radial direction. The parameter a in (j3)) governs the evolution of the extra dimension and serves as 
a measure of the anisotropy. Namely, for a = ±1 the fluid has isotropic pressure and the equation 
of state is that of radiation p = 3p, whereas for any other a the fluid has anisotropic pressures 
obeying p = p + 2p^. 

An observer moving in the radial direction with velocity v given by (|T8|) will interpret the matter 
distribution as an isotropic fluid satisfying p = 3p, with a radial heat flux directed outward, for any 
value of a in the regions < a < 1 and -^3 < a < -1. The energy density and the temperature 
are related by the Stefan-Boltzmann law (j2Hl). We note that the introduction of energy flux is 
crucial here to account for Tqi ^ 0. This is different from other solutions in the literature where 
Toi 7^ can be explained in terms of radial motion only. 

The effective energy-momentum tensor satisfies the appropriate energy conditions but diverges 
at the hypersurfaces rt = 0. From a 5D viewpoint, the latter is associated with the singularity of 
the Kretschmann scalar. But, from a AD viewpoint it is related to the homothetic nature of the 
distribution, which is singled out from a complicated set of initial conditions. 

An important feature here is that the matter density and pressure are spatially nonuniform. 
This suggest that our solution can be applied to describe the evolution of a spherical inhomogeneity 
in an early universe. 

We would like to finish with the following remarks. Firstly, that although the 4Z}-solution 
was obtained in the context of STM in 5D, by virtue of Campbell's theorem it is also a solution of 
the AD Einstein field equations with energy- momentum Secondly, the same matter distribution 
is obtained if the metric Q is interpreted in the context of brane theory in 5D ^E] . Finally, the 
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metric © exhibits a number of interesting properties, which make it different from other solutions 
in gravity. We will discuss this in more detail elsewhere. 
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